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Abstract. Global mapping of terrestrial vegetation fluores- 
cence from space has recently been accomplished with high 
spectral resolution (v/Av > 35 000) measurements from the 
Japanese Greenhouse gases Observing SATellite (GOSAT). 
These data are of interest because they can potentially pro- 
vide global information on the functional status of vegeta- 
tion including light-use efficiency and global primary pro- 
ductivity that can be used for global carbon cycle model- 
ing. Quantifying the impact of fluorescence on the O 2 -A 
band is important as this band is used for photon pathlength 
characterization in cloud- and aerosol-contaminated pixels 
for trace-gas retrievals including CO 2 . Here, we examine 
whether fluorescence information can be derived from space 
using potentially lower-cost hyperspectral instrumentation, 
i.e., more than an order of magnitude less spectral resolu- 
tion (v/Av~1600) than GOSAT, with a relatively simple 
algorithm. We discuss laboratory measurements of fluores- 
cence near one of the few wide and deep solar Fraunhofer 
lines in the long-wave tail of the fluorescence emission re- 
gion, the calcium (Ca) II line at 866 nm that is observable 
with a spectral resolution of ~0.5 nm. The filling-in of the 
Ca II line due to additive signals from various atmospheric 
and terrestrial effects, including fluorescence, is simulated. 


We then examine filling-in of this line using the SCanning 
Imaging Absorption spectroMeter for Atmospheric CHar- 
tographY (SCIAMACHY) satellite instrument. In order to 
interpret the satellite measurements, we developed a gen- 
eral approach to correct for various instrumental artifacts 
that produce false filling-in of solar lines in satellite mea- 
surements. The approach is applied to SCIAMACHY at the 
866 nm Ca II line and to GOSAT at 758 and 770 nm on the 
shoulders of the O 2 -A feature where there are several strong 
solar Fraunhofer lines that are filled in primarily by vegeta- 
tion fluorescence. Finally, we compare temporal and spatial 
variations of SCIAMACHY additive signals with those of 
GOSAT and the Enhanced Vegetation Index (EVI) from the 
MODerate-resolution Imaging Spectroradiometer (MODIS). 
Although the derived additive signals from SCIAMACHY 
are extremely weak at 866 nm, their spatial and temporal 
variations are consistent with chlorophyll a fluorescence or 
another vegetation-related source. We also show that filling- 
in occurs at 866 nm over some barren areas, possibly origi- 
nating from luminescent minerals in rock and soil. 
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1 Introduction 

Fluorescence from terrestrial vegetation has been measured 
in the laboratory and with ground- and aircraft-based instru- 
ments for several decades (see e.g., the review of Meroni et 
al., 2009, and the many references therein). More recently, it 
has been demonstrated that fluorescence information can also 
be derived from satellite passive sensors measuring backscat- 
tered sunlight at wavelengths near the O 2 -A band absorption 
feature. It is important to quantify the impact of fluores- 
cence on atmospheric absorption bands, such as the O 2 -A 
band, that are used for cloud- and aerosol-characterization 
for other trace-gas retrievals, including CO 2 (e.g., Franken- 
berg et al., 2011a). These fluorescence measurements, par- 
ticularly with global coverage from satellites, are of inter- 
est in their own right because the signal originates from the 
core complexes of the photosynthetic machinery where en- 
ergy conversion of absorbed photosynthetically active radia- 
tion (PAR) occurs. Chlorophyll fluorescence measurements 
therefore may provide a means to estimate global instanta- 
neous vegetation carbon-related processes, such as those im- 
portant to precision farming, forest management, and assess- 
ment of the terrestrial carbon budget including gross primary 
productivity (GPP) (e.g., Campbell et al., 2008; Damm et al., 
2010; Frankenberg et al., 2011b). 

Fluorescence in vegetation occurs because a portion of the 
sunlight absorbed by chlorophyll cannot be used for carbon 
fixation. While most of this unused energy is given off as 
heat, a small fraction is re-emitted as fluorescence at longer 
wavelengths (lower energy) with respect to the excitation 
wavelength. The far-red and red fluorescence generally oc- 
curs as a convolution of broad band emission with two peaks 
at 685 and 740 nm, respectively, as shown in Fig. 1 (e.g., 
Meroni et al., 2009; Corp et ah, 2003, 2006). Many stud- 
ies have shown that in high light conditions (i.e., in the late 
morning and early afternoon when many satellite measure- 
ments are made) and when plants are under stress, fluores- 
cence is correlated with photosynthesis and light use effi- 
ciency (LUE) (e.g., Flexas et ah, 2002; Louis et ah, 2005; 
Meroni et ah, 2008; van der Tol et ah, 2009; Zarco-Tejada 
et ah, 2009; Daumard et ah, 2010; Amoros-Lopez et ah, 
2008). Other research has indicated that chlorophyll fluo- 
rescence supplies information content about photosynthetic 
function that is complementary to reflectance -based spectral 
vegetation indices including the Normalized Difference Veg- 
etation Index (NDVI) and the Enhanced Vegetation Index 
(EVI), linked to chlorophyll content, and the Photochemical 
Reflectance Index (PRI), related to changes in xanthophyll 
cycle pigments (Meroni and Colombo, 2006; Middleton et 
ah, 2008, 2009; Rascher et ah, 2009; Meroni et ah, 2008; 
Daumard et ah, 2010; Guanter et ah, 2007; Zarco-Tejada 
et ah, 2009; Joiner et ah, 2011; Frankenberg et ah, 2011b; 
Guanter et ah, 2012). 

Space-based measurement of solar-induced chlorophyll 
fluorescence is challenging, because its signal (typically 


Fig. 1. Measured fluorescence spectra of excised com leaves (Zea 
Mays L.) in the laboratory, with simulated solar illumination at full 
mid-day sun intensity as measured below a long-pass filter blocking 
wavelengths longer than 665 nm, showing maximum fluorescence 
(Fmax) after 5 min. of dark adaptation, steady state fluorescence 
(F s , achieved in less than 0.5 min after illumination), and a sample 
solar spectrum from Thuillier et al. (2003, 2004). 


1-5 % at red and far-red wavelengths) is small as compared 
with the much larger reflectance signal. Ground- and aircraft- 
based approaches have made use of the dark and spectrally- 
wide O 2 -A (~760nm) and O 2 -B (~690 nm) atmospheric 
features to detect the weak fluorescence signal (see e.g., 
Meroni et ah, 2009). Deep solar Fraunhofer lines have also 
been used to detect fluorescence from vegetation using air- 
borne and ground-based sensors. Initial studies (e.g., Plas- 
cyk and Gabriel, 1975; McFarlane etah, 1980; Watson, 1981; 
Sioris et ah, 2003) focused on lines shortward of the red edge 
(~700nm) where the terrestrial surface reflectivity is rela- 
tively low. More recently. Joiner et ah (2011), Frankenberg 
et ah (2011b), and Guanter et ah (2012) focused on longer- 
wavelength solar Fraunhofer lines that can be observed with 
space -based instruments such as the currently operational 
GOSAT. They showed that fluorescence can be detected us- 
ing Fraunhofer lines off the far-red chlorophyll a fluores- 
cence peak even when the surface is relatively bright. 

Here we build on that work by examining the possibil- 
ity of making fluorescence measurements from space with 
lower spectral resolution instrumentation as compared with 
the GOSAT interferometer or similar instruments. If this can 
be demonstrated, then low cost instrumentation can poten- 
tially be used for space- and air-borne measurements with a 
relatively simple retrieval algorithm, i.e., one that does not 
require complex atmospheric correction. Of the wide and 
deep solar lines shown in Fig. 1 , we focus on the 866 nm Ca II 
solar Fraunhofer line. This is the only line with negligible 
contamination from atmospheric water vapor capable of de- 
tecting fluorescence at spectral resolutions of about 0.5 nm. 
We first examine laboratory measurements of fluorescence at 
the long-wavelength tail of the 740 nm feature near 866 nm. 
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We then simulate the potential filling-in of this line from at- 
mospheric effects, namely vibrational and rotational Raman 
scattering, as well as from fluorescence and high temperature 
emitters such as fires and volcanoes. Finally, we make use of 
unique spectral radiance measurements from the SCanning 
Imaging Absorption spectroMeter for Atmospheric CHar- 
tographY (SCIAMACHY) instrument aboard the European 
Space Agency (ESA) Environmental satellite (EnviSat) en- 
compassing the Ca II line at 866 nm. SCIAMACHY obser- 
vations are used to retrieve an additive signal that fills in the 
Ca II line. This is clearly not an optimal wavelength region 
for measuring fluorescence owing to the extreme weakness 
of the signal; instrumental artifacts could easily overwhelm 
a fluorescence signal. However, although not optimized for 
detecting fluorescence, SCIAMACHY can be used for an ex- 
ploratory study. 

An integral part of this work is to introduce a general ap- 
proach to correct for a number of different instrumental arti- 
facts that produce false filling-in signals and thus bias flu- 
orescence retrievals. The methodology is applied to both 
SCIAMACHY and GOS AT data and is generalizable to other 
instruments. This methodology improves the accuracy of the 
results reported in the initial GOSAT study of Joiner et al. 
(2011) and enables the use of more data as compared with 
the approach of Frankenberg et al. (2011b), thereby poten- 
tially increasing the fidelity of GOSAT data. 

The paper is organized as follows: we discuss laboratory 
fluorescence measurements at the long-wavelength tail of the 
chlorophyll a fluorescence spectrum in Sect. 2. Section 3 de- 
scribes the satellite observations used here. Simulations of 
the satellite-observed Ca II line-filling at 866 nm from at- 
mospheric and terrestrial effects including fluorescence are 
given in Sect. 4. In this section, we also examine the impact 
of clouds and aerosol on fluorescence retrievals. The details 
of our retrieval approach, including methodology to correct 
for instrumental artifacts, are provided in Sect. 5. We show 
maps of monthly-mean derived additive signals from 866 nm 
Ca II line -filling and a seasonal climatology using more than 
8 yr of SCIAMACHY data in Sect. 6. In this section, we 
also compare and contrast temporal and spatial variations 
of the additive signals (retrieved from both SCIAMACHY 
and GOSAT) with the MODIS EVI. Conclusions are given 
in Sect. 7. 


2 Laboratory measurements of fluorescence at 
wavelengths >850 nm 

As the derived additive signals at 866 nm shown in Sect. 6 
appear to be consistent with a vegetation source, namely flu- 
orescence, here we discuss laboratory measurements of veg- 
etation fluorescence made at similar wavelengths. Labora- 
tory instruments using light-induced pulses have focused on 
the red and far-red fluorescence emission peaks and not the 
far-red/NIR emission tail in diagnostic measurements. This 


is true for laser-induced lab and field measurements as well. 
Detectors are typically chosen to be affordable and are there- 
fore usually limited in spectral range. Consequently, very 
little information is available on the NIR tail region beyond 
about 800 nm. 

Figure 1 displays fluorescence measurements of corn 
leaves acquired in our laboratory using polychromatic exci- 
tation at wavelengths shorter than 665 nm (Oriel 91 160A and 
Oriel 81080) after a five minute dark adaptation of the sam- 
ples. Chlorophyll fluorescence levels decrease significantly 
from the initial dark-adapted maximum ( F max ) to the steady 
state (F s ) in 15-20 s, after which they remain relatively con- 
stant over time. Here, we show that F max and F s levels in the 
red and far-red peaks vary significantly, but are much lower 
in value and do not differ significantly in the emission tails 
(e.g., at 656 and 854-866 nm). At 866 nm, the measured sig- 
nal is of the order of 0. 1-0.2 mW m~~ nm -1 sr“ 1 . 

Amoros-Lopez et al. (2008) similarly report spectral mea- 
surements from 650-850 nm showing a small amount of fluo- 
rescence (~0.1 mW m -2 nm -1 sr -1 ) at 850 nm when illumi- 
nated by sunlight through a filter that restricts excitation to 
wavelengths shorter than about 600 nm (i.e., not full ambient 
intensity at leaf level). They also measured a diurnal cycle 
of fluorescence for wavelengths up to 800 nm using artificial 
illumination in a setup similar to that described above. Their 
results show very little difference in fluorescence measured 
at 10:00 (close to the SCIAMACHY local overpass time) and 
12:00 (closer to the GOSAT local overpass time) for wave- 
lengths between 770 and 800 nm. Fluorescence at 14:00 was 
slightly lower (by ~15%) as compared with noon. Other 
measurements (e.g., Chappelle and Williams, 1987; Rosema 
et al., 1991; Edner et al., 1994; Saito et al., 1998) suggest that 
the long-wave tail of the chlorophyll a fluorescence emis- 
sions and perhaps fluorescence from other compounds such 
as vitamin Ki (Chappelle et al., 1984) may produce a mea- 
surable signal at wavelengths beyond 800 nm. 

We did not, however, detect a resolvable fluorescence sig- 
nal with respect to instrument signal-to-noise at 866 nm us- 
ing a Fluorolog® 3 spectrofluorometer with monochromatic 
excitation. Excitation wavelengths ranging from the ultravi- 
olet through the visible were used. One possible explana- 
tion is that fluorescence at the long wavelength tail of the 
fluorescence spectrum may be higher when induced by fiill- 
spectrum sunlight as compared with a monochromatic source 
or in a partially-illuminated scenario. Measurements have 
shown that the excitation efficiency of chlorophyll fluores- 
cence is high in the blue and ultraviolet (UV) with a mini- 
mum around 550 nm and another peak of high efficiency be- 
yond about 650 nm (Chappelle and Williams, 1987; Rosema 
et al., 1991; Middleton et al., 2008). 
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3 Satellite observations 

3.1 SCIAMACHY 

SCIAMACHY is a grating spectrometer that measures trans- 
mitted, reflected, and scattered sunlight in both limb- and 
nadir-viewing geometries from ultraviolet to near-infrared 
wavelengths (212-2386 nm) in eight separate channels. It 
was launched on the EnviSat platform in February 2002 
into a sunsynchronous orbit with a descending node equa- 
tor crossing time near 10:00 LT. The instrument can make 
measurements in a number of different modes (e.g., Licht- 
enberg et ah, 2006). In addition to Earth views in the limb 
and nadir, it can operate in a solar or lunar occultation mode. 
It also measures the solar irradiance. In this work, we pri- 
marily use channel 5 that covers wavelengths between 773 
and 1063 nm at a spectral resolution of 0.54 nm. The signal- 
to-noise ratio (SNR) for nadir measurements in this channel 
is between 1000 and 10 000 (Gottwalk et al., 2006) and the 
ground footprint size is approximately 30 km in the along 
track direction and about 60 km in the cross-track direc- 
tion for latitudes between 60° N and 60° S. We also exam- 
ined spectra in channel 4 that covers 595-812 nm at 0.48 nm 
spectral resolution, but do not include any results from that 
channel in this work. 

We processed raw SCIAMACHY level IB data with the 
SciaLlc command-line tool software package (DLR, 2006). 
Using this package, we applied all available corrections and 
calibrations. Although the corrections do an excellent job of 
removing most artifacts that produce false filling-in of solar 
Fraunhofer lines, we found it necessary to further account for 
the remaining effects as discussed in Sect. 5. In Appendix A, 
we give a brief overview of the effects. More details can be 
found in Lichtenberg et al. (2006). 

3.2 GOSAT 

GOSAT is a satellite mission designed to monitor the global 
distribution of the greenhouse gases CO 2 and CH 4 (Yokota 
et al., 2009). Jointly developed by the Japanese Ministry 
of the Environment (MOE), the National Institute for Envi- 
ronmental Studies (NIES), and the Japanese Aerospace Ex- 
ploration Agency (JAXA), GOSAT was launched on 23 Jan- 
uary 2009 into a sun-synchronous orbit. It has a descending 
node equatorial crossing time near 13:00 LT and a 3 day re- 
peat cycle. GOSAT has two instrument packages: (1) the 
Thermal And Near-infrared Sensor for carbon Observation- 
Fourier Transform Spectrometer (TANSO-FTS) (Kuze et al., 
2009) measures backscattered solar radiation in three short- 
wave infrared (SWIR) regions, referred to as “bands”, cen- 
tered at 0.76, 1.6, and 2.0 pm with a nadir ground footprint 
of 10.5 km diameter and (2) the Cloud and Aerosol Imager 
(CAI) contains four bands at 0.38, 0.67, 0.87, and 1.6 pm 
and footprints between 0.5 and 1.5 km. 


Joiner et al. (2011), Frankenberg et al. (2011b), and Guan- 
ter et al. ( 2012 ) retrieved chlorophyll fluorescence using 
TANSO-FTS band 1. This band extends from approximately 
758 to 775 nm with a resolving power (v/Av) of >35 000. 
The SNR for this band is >300 for a typical scene radiance. 

There are also instrumental effects that can distort GOSAT 
spectra to produce a false filling-in of Fraunhofer lines. 
These include non-linearity from the analog-to-digital con- 
verter (ADC) (Kuze et al., 2011b) and the analog filter cir- 
cuit in band 1 that can lead to effects such as a zero-level off- 
set (Frankenberg et al., 2011b). In addition, calibration drift 
(Kuze et al., 2011a) may produce errors and drift in the ab- 
solute values of the derived additive signal. The latest avail- 
able version (vl30) of the level lb data addresses some, but 
not all of these problems. Future versions will provide more 
corrections for these non-linearities. Only a limited amount 
of vl30 data is currently available. Therefore, here we use 
a previous version of level lb data (v050-vl00) with an ap- 
proach to account for instrumental artifacts that produce false 
filling-in as described below. 

3.3 MODIS Vegetation indices 

Below, we compare satellite-derived fluorescence with the 
Enhanced Vegetation Index (EVI) (Huete et al., 2002), a 
popular vegetation reflectance -based index that indicates rel- 
ative greenness and is used to infer photosynthetic func- 
tion. The EVI, derived from the MODIS instrument, can be 
found in the “Vegetation Indices 16-Day L3 Global 1km” 
(MOD13A2) data set. We use data from the Aqua MODIS, 
acquired in the early afternoon, exclusively here for direct 
comparisons with 2009 GOSAT results and averaged over 
the period August 2002 through June 2011 for climatologi- 
cal comparisons with SCIAMACHY. The Aqua satellite has 
an ascending node equator crossing near 13:30 LT. Although 
these data have been cloud-filtered and corrected for aerosol- 
corruption effects, we followed the methodology of Xu et al. 
( 2011 ) to further remove residual atmospheric corruption in 
order to produce the best possible signals from vegetation. 

4 Simulated filling-in at 866 nm 

Here, we simulate the effects of additive signals such as flu- 
orescence as well as rotational- and vibration-Raman scat- 
tering on space-based observations of the Ca II line near 
866 nm. We use the so-called KPNO2010 high spectral res- 
olution solar irradiance reference spectrum (Chance and Ku- 
rucz, 2010 ) for these simulations and convolve the spectrum 
with triangular instrument line shape functions having differ- 
ent values of the full-width half-maximum (FWHM) includ- 
ing one similar to SCIAMACHY. 
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Fig. 2. Simulated spectral filling-in around the Ca II solar line from 
an additive signal of 0.2 mW m -2 nm" 1 sr - 1 , SZA = 40° nadir, and 
surface reflectivity = 0.3 (RRS not included) for different instrument 
resolutions (triangular slit function) with full-width half-maximum 
(FWHM) as indicated including that of SCIAMACHY (SCIA). 


4.1 Fluorescence and other additive signals 

Figure 2 shows simulated filling-in of the 866 nm Ca II 
solar line for a spectrally constant additive signal of 
0.2 mWm -2 sr -1 nm -1 , for several different spectral reso- 
lutions including the limit at the native 0.01 nm resolution of 
the Chance and Kurucz (2010) spectrum (denoted monochro- 
matic in the figure). Because Rayleigh optical thickness is 
low at these wavelengths, we did not include atmospheric 
scattering in these simulations. The surface is treated as 
Lambertian with an equivalent reflectivity. The filling-in at 
SCIAMACHY resolution is small (<1 %). However, modest 
increases in spectral resolution will increase the filling-in. 

Figure 3 shows the surface reflectivity dependence of the 
simulated filling-in of the 866 nm Ca II line for the same sig- 
nal (0.2 mWm -2 nm -1 sr -1 ) and spectral resolutions as in 
Fig. 2. The surface reflectivity dependence is significant at 
low surface reflectivities. However, the dependence flattens 
out at values typical of vegetated land (>~0.25). 

Note that the above results hold for any type of additive 
signal including fluorescence from sources in vegetation or 
minerals. Another potential source of additive signal is radi- 
ation emitted from fires or volcanoes. Using the Planck func- 
tion, we estimate that the radiance from a source at 1000 K 
will produce a signal of the order of 15 mW m~ 2 nm -1 sr -1 
at 866 nm and ~3 mW m“ 2 nm -1 sr~' at 770 nm. If such 
a source filled a SCIAMACHY or GOSAT field-of-view, it 
would certainly be detectable. However, flaming portions 
of fires and sources of hot gas and lava from volcanoes are 
most likely to fill only a small fraction of the relatively large 
GOSAT or SCIAMACHY footprints. Never-the-less, fires 
and volcanoes should be considered as potential sources of 
solar line filling at far-red and NIR wavelengths. Reducing 



Surface albedo 


Fig. 3. Filling-in due to an additive signal of 0.2 mW m -2 nm -1 
sr - , SZA = 20°, and nadir as a function of surface reflectivity for 
different instrument spectral resolutions. 


temperatures to <~650 K renders the filling-in signals unde- 
tectable at the wavelengths and resolutions of SCIAMACHY 
and GOSAT. 

4.2 Rotational-Raman scattering 

We computed the filling-in of the 866 nm Ca II line due 
to rotational-Raman scattering (RRS) using the scalar Lin- 
earized Discrete Ordinate Radiative Transfer code (LIDORT- 
RRS) of Spurr et al. (2008). Figure 4 shows the filling-in as 
in Fig. 2 but from rotational-Raman scattering (RRS) only. 
Although the amount of atmospheric scattering is relatively 
small at these wavelengths, the RRS filling-in is apparent ow- 
ing to the large depth of this solar line as similarly shown by 
Joiner et al. (2011) for the K I line at 770 nm. The spec- 
tral signature is very similar to that of line filling from a 
spectrally flat additive signal. The filling-in owing to RRS 
is about a factor of 6 less than that due to an additive signal 
of 0.2 mW m -2 nm -1 sr -1 at SCIAMACHY resolution. This 
factor increases somewhat at higher spectral resolutions. 

Figure 5 shows that the RRS filling-in dependence on 
surface reflectivity is relatively flat for reflectivities >~0.3. 
There is a small but significant dependence on both surface 
pressure and solar zenith angle (SZA). Figure 6 shows the 
SZA dependence in more detail as well as the dependence 
on satellite view zenith angle (VZA) and surface pressure. 
The SZA dependence is relatively small for SZA < 60°, but 
increases rapidly for greater angles. The VZA dependence is 
small for typical angles viewed by SCIAMACHY. 

4.3 Vibrational-Raman scattering 

Vibrational Raman scattering (VRS) is commonly neglected 
in UV and visible retrieval algorithms as it is small compared 
with rotational-Raman scattering (RRS) (Burrows et al., 
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Fig. 4. Simulated filling-in spectra due to rotational-Raman scat- 
tering (RRS) only at SZA = 40°, nadir, and surface reflectivity = 0.3 
for different instrument spectral resolutions. 
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Fig. 5. Filling- in of Ca II line at 866 nm due to rotational-Raman 
scattering as a function of surface reflectivity for two solar zenith 
angles (SZA) and two surface pressures (P s ) showing simulated dif- 
ference between low clouds over ocean (800 hPa) and a low lying 
land area (1000 hPa). 

1996). However, the fractional contribution of VRS is larger 
in the NIR owing to higher values of solar radiation at exci- 
tation wavelengths as compared with emission wavelengths. 
The VRS frequency shifts between excitation and emission 
are An ~ 2331 and 1555 cm -1 for N 2 and O 2 , respectively, 
assuming pure vibrational transitions of the molecules (Heise 
and Schrotter, 1995). We assessed the VRS contribution 
to the filling-in of the 866 nm Ca II Fraunhofer line using 
the single scattering approximation of the radiative transfer 
equation. Along with photons only scattered in the atmo- 
sphere, we also accounted for photons first reflected from the 
surface and then scattered in the atmosphere. 

Absolute differential Raman cross sections of the vibra- 
tional line at 488 mn at 90° measured by different authors 
for N 2 vary from 3.3 x 10 -31 to 5.6 x 10 -31 cm 2 str-mol -1 



SZA (deg.) 

Fig. 6. Filling-in of Ca II line at 866 nm due to rotational-Raman 
scattering as a function of solar zenith angle (SZA) for 2 view zenith 
angles (VZA) typical of minimum (nadir, 0°) and maximum (40°) 
for SCIAMACHY for a relative azimuth angle (AZ) of 90°. 

(Fenner et al., 1973). Here, we use a value of 
5.5 x 10 - ’ 1 cm 2 str-mol -1 that is consistent with data from 
Heise and Schrotter (1995). The ratio of O 2 to N 2 cross sec- 
tions at 488 nm is assumed to be 1.3 (Fenner et al., 1973). 
We adopt a general form of phase functions for both the un- 
shifted Cabannes line and wavelength shifted VRS lines as 
suggested by Danichkin et al. (1982). For the central Ca- 
bannes line in air, the depolarization ratio for unpolarized 
light is assumed equal to 0.00734 (Chance and Spurr, 1997). 
For VRS lines, we use values recommended by Danichkin et 
al. (1982), i.e., 0.043 and 0.099 forN 2 andC> 2 , respectively. 

Figure 7 shows the VRS filling-in of the 866 nm Ca II 
Fraunhofer line as a function of surface reflectivity. The 
maximum value of about 0.45 % occurs for a non-reflecting 
surface. VRS filling-in rapidly decreases with increasing sur- 
face reflectivity. For typical values of surface reflectivity over 
land (>0.2), the filling-in only weakly depends on surface re- 
flectivity and does not exceed 0.02 % for SZA < 60°. There- 
fore, we neglect VRS in the remainder of this study. 


5 Retrieval methodology 
5.1 Fitting windows 

The use of solar Fraunhofer lines substantially simplifies 
a space -based retrieval of F s as compared with the O 2 - 
A band as described by Joiner et al. (2011) and Franken- 
berg et al. (2011a); atmospheric correction is unnecessary 
as long as observations are made at wavelengths not af- 
fected by atmospheric absorption. Joiner et al. (2011) used 
the potassium (K) I line near 770 nm to derive an addi- 
tive signal using GOSAT’s TANSO-FTS. This line falls be- 
tween several weak O 2 lines and can be observed in relative 
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Surface albedo 

Fig. 7. Filling-in of Ca II line at 866 nm due to vibrational-Raman 
scattering (VRS) as a function of surface reflectivity for 2 solar 
zenith angles (SZA). 


isolation with GOSAT’s high spectral resolution. Here, we 
use the same GOSAT fitting window as Joiner et al. (2011) 
(769.90-770.25 nm). We use a second fitting window be- 
tween 758.45 and 758.85 nm. Our fitting windows are simi- 
lar to, though more narrow than, those used by Frankenberg 
et al. (2011b) and Guanter et al. (2012). While the use of 
broader fitting windows that include more Fraunhofer lines 
may potentially improve the signal-to-noise ratio of the re- 
trievals, it may also require the retrieval of additional param- 
eters such as those that account for the wavelength depen- 
dence of surface reflectivity and fluorescence. We have cho- 
sen here a conservative approach that utilizes narrow win- 
dows that contain the strongest Fraunhofer lines, although 
the approach can be extended to wider fitting windows. Fur- 
ther optimization of fitting windows is beyond the scope of 
this work but may be the focus of a future study. 

Of the solar lines shown in Fig. 1 (656, 850, 854, 866 nm), 
only the 866 nm Ca II line is relatively free of contamina- 
tion by atmospheric water vapor absorption. We attempted 
to retrieve additive signals using all of these lines including 
the Ha line near 656 nm in SCIAMACHY channel 4. We 
found that it was difficult to disentangle the filling-in signal 
from water vapor absorption. Therefore, in the remainder of 
this work, we focus solely on deriving an additive signal us- 
ing the 866 nm Ca II solar line from SCIAMACHY with the 
spectral band 863.5-868.5 nm. 

5.2 General approach 

We use a variation of the approach implemented by Joiner 
et al. (2011) to jointly retrieve a reflectivity-related param- 
eter K and an additive signal. We denote the generic addi- 
tive signal henceforth as “F”. We use the above-mentioned 
spectral fitting windows with the following simplified model 
for the observed Earth spectral radiance I(X) that assumes 


negligible atmospheric absorption and scattering (Joiner et 
al., 2011): 

7(1) = K E(X)* + F, (1) 

where E(X)* is a reference spectrum ideally containing no 
filling-in from the source of interest. We note that if the ref- 
erence spectrum contains filling-in from other sources, the 
retrieved value will reflect a residual additive signal from 
sources not present in the reference. 

As in Joiner et al. (2011), we assume no wavelength de- 
pendence of F over our relatively small fitting windows. We 
fit both constant and linear in wavelength coefficients for K 
(i.e., K = Kq + K\ [X — ko]) for the larger SCIAMACHY fit- 
ting window, but only a constant for the smaller GOSAT win- 
dows. The spectral dependence of atmospheric, cloud, and 
aerosol scattering is essentially constant over the relatively 
small narrow fitting windows used, so the neglect of it will 
have an insignificant effect on the retrieval of F. Similarly, 
calibration of the reference spectrum E(X)* does not affect 
the derived value of F, since any offset between the reference 
and observed spectra is absorbed in the parameter K. An er- 
ror in absolute calibration of 7(k), however, will produce an 
error in F. 

Similar to the approach of Joiner et al. (2011), we use a 
standard weighted least squares fitting procedure. All wave- 
lengths are weighted equally except for a few wavelengths 
within the GOSAT fitting window that are not used because 
they are affected by O 2 absorption. 

5.3 Use of reference spectra 

The main difference between our approach and that of Joiner 
et al. (2011) is that we use spectral radiance measurements 
made over the cloudy ocean as a reference rather than a mea- 
sured or computed solar irradiance spectrum. No further cor- 
rections to observed spectra are made. Our retrieval approach 
therefore does not require a reference solar spectra. In addi- 
tion, it is relatively insensitive to errors in the instrument line 
shape function (ILSF), to the extent that the ILSF is suffi- 
ciently stable. 

The main criterion for reference spectra used to derive F 
from a terrestrial source is that the reference spectra should 
represent the spectral shape of a solar Fraunhofer line with- 
out filling-in from the terrestrial source. There are several 
advantages of using cloudy Earth radiance spectra as a ref- 
erence to derive a terrestrial additive signal as opposed to 
solar irradiance spectra (measured or computed). Firstly, the 
use of Earth spectra minimizes the so-called undersampling 
problem and eliminates errors due to inaccuracy of the as- 
sumed instrument line shape function. The undersampling 
problem and a proposed correction method have been dis- 
cussed at length by Chance et al. (2005). However, we found 
that the proposed correction method did not achieve the de- 
sired accuracy owing to uncertainties in the instrument line 
shape function. Because the wavelength shifts are very small 
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between ocean reference spectra and the land spectra of in- 
terest, problems due to undersampling are essentially elimi- 
nated. Corrections used by Joiner et al. (2011), that account 
for under-sampling when referencing to observed solar spec- 
tra slightly shifted in wavelength as compared with Earth 
spectra, are not necessary. 

Secondly, referencing to ocean spectra with similar con- 
tinuum radiance values allows us to significantly reduce the 
errors from stray light, dark current, zero-level offset, and 
other instrumental effects that can produce a false filling-in. 
These effects should be present in the ocean reference spectra 
at approximately the same level as in land spectra. Therefore, 
the effects should be accounted for when deriving a terrestrial 
additive signal using Eq. (1). 

Finally, using ocean spectra at similar radiance levels and 
SZAs as a reference will implicitly account for some of 
the filling-in due to RRS. RRS depends upon surface re- 
flectivity as well as the mean pressure of the scattering sur- 
face. Flere, we are using cloudy ocean data at similar SZAs 
and continuum radiance values to land spectra. Studies of 
the so-called optical centroid pressure (approximately, the 
reflectance-weighted mean pressure) of clouds over ocean at 
moderate reflectivities typical of vegetated land show that 
the clouds occur mainly between about 775 and 875 hPa 
(Joiner et al., 2012). Deep convective and frontal clouds 
that have lower optical centroid pressures typically occur 
only at higher Lambertian-equivalent reflectivities (>~0.6). 
Therefore, the main difference in RRS filling-in between the 
cloudy ocean and land spectra is due to the different effec- 
tive surface pressures. As shown in Fig. 5, this difference 
is small, but not negligible for the expected pressure differ- 
ences, sometimes leading to negative retrieved values of F. 
We must bear this difference in mind when examining our 
results. For regions with surface pressures higher than about 
900 hPa, there may be unaccounted - for filling - in due to 
RRS. However, for high-altitude land regions (<775 hPa) the 
reference ocean spectra may overestimate RRS filling-in. 

Frankenberg et al. (2011b) used a different approach to ac- 
count for instrumental artifacts; they utilized measurements 
over Antarctica to develop a parameterization for the false 
fluorescence produced by the zero-level offset in GOSAT 
data. They showed generally increasing amounts of false 
fluorescence with radiance levels with some wiggles and a 
downturn at the highest radiances in the S polarization. We 
found similar features using cloudy ocean data. With our 
approach, we are able to cover the full range of observed ra- 
diance values and see a downturn in the false fluorescence 
or filling-in at the highest radiance levels also in the P po- 
larization. Our utilization of cloudy ocean data does not ne- 
cessitate the use of high SZA data to cover the full range of 
observed continuum radiance values; filling-in from RRS in- 
creases rapidly at high SZAs as shown above, although the 
imprint of RRS is not obvious in the results of Frankenberg 
et al. (2011b). 


Frankenberg et al. (2011b) discarded GOSAT P polariza- 
tion data because the errors in that polarization were time- 
dependent. We likewise find such time-dependent errors in 
both GOSAT (both polarizations) and SCIAMACHY data as 
shown in Appendix B. Because cloudy ocean data are avail- 
able year round, we are able to account for these errors and 
use all data including GOSAT P polarization, thereby signif- 
icantly reducing the effects of instrumental noise. 

In addition, cloudy ocean data span all latitudes. There- 
fore, we are able to account for instrumental effects that 
vary systematically with orbital phase; our approach is able 
to detect and correct for subtle but significant latitudinally- 
dependent artifacts in both GOSAT and SCIAMACHY data. 
We find that the latitudinal-dependence of these artifacts also 
varies with season. Further details can be found in Ap- 
pendix B. Guanter et al. (2012) similarly use data over cloudy 
ocean to account for the GOSAT zero-level offset, but appar- 
ently do not account for latitudinal dependences as we have 
here. 

Here, we binned potential reference spectra according to 
several criteria discussed below. For each bin, we averaged 
multiple spectra to produce a composite reference spectrum. 
To reduce any filling-in signal from the ocean, e.g., from 
ocean Raman scattering (Vasilkov et al., 2002) or fluores- 
cence, we only use observations with reflectances >0.2 for 
SCIAMACHY. For GOSAT we used only observations with 
a cloud fraction (derived from the CAI) of >20%. 

5.3.1 SCIAMACHY binning 

For SCIAMACHY, we binned spectra separately for each of 
the two primary pixel exposure times and for each scan posi- 
tion within a nadir block according to the continuum radiance 
value and the wavelength shift (with respect to a fixed solar 
spectrum). We found that the observed filling-in over ocean 
sometimes changes abruptly and that it also drifts slowly over 
time. We therefore compute a set of reference spectra for 
each day of observations. If the standard deviation of the 
filling-in over ocean for a given day is high, we eliminate 
that day from our sample. 

We use a bin size of 0.1 for pcos(SZA) (proportional to 
normalized radiance) where p is the reflectance. We found 
that the wavelength shifts varied systematically with latitude, 
so binning by both wavelength shift and reflectance implic- 
itly accounts for filling-in features such as RRS that are func- 
tion of solar zenith angle as well as instrumental artifacts that 
vary systematically with orbital phase. Because the wave- 
length shifts vary with time, we compute the daily mean and 
standard deviation of wavelength shift and use 8 bins with in- 
tervals of half the standard deviation of the wavelength shift 
(typically ~0.007 nm). 


Atmos. Meas. Tech., 5, 809-829, 2012 


www.atmos-meas-tech.net/5/809/2012/ 


J. Joiner et al.: Filling-in of near-IR solar lines 


817 


5.3.2 GOSAT binning and additional corrections 

For GOSAT, we used bin sizes of 0.05 for normalized ra- 
diance 20° for latitude, and 0.001 nm for wavelength shift. 
We use only data processed with the H gain range as there 
are not enough data from the M range to perform a reliable 
correction. Therefore, data over the Sahara and in parts of 
central Australia have not been processed. For GOSAT, we 
retrieve F separately for the S and P polarizations (separate 
corrections applied) and average them. We note small differ- 
ences in absolute values obtained with the two polarizations. 
Finally, we apply a correction factor to the radiances to ac- 
count for instrument degradation as suggested by Kuze et al. 
(2011a). 

5.4 Cloud filtering 


Because clouds do not alter the spectral structure used to re- 
trieve F, we can tolerate small to moderate amounts of cloud 
within our satellite footprints; the main effect of clouds will 
be to shield a portion of the terrestrial additive signal from 
satellite view. We remove significantly cloud-contaminated 
data from our samples by two different methods. For SCIA- 
MACF1Y, we use the concept of the cloud radiance fraction 
(CRF). CRF is defined as the fraction of observed radiance 
that is scattered by cloud particles. The CRF is a convenient 
quantity for cloud screening, because it can be computed di- 
rectly for large fields-of-view and collocation with a subpixel 
imager is not necessary. 

The CRF is given by 


CRF = / g 


fcld 

fobs 


(2) 


MODIS (MOD43B3) at 865 nm (Lucht et al., 2000). We 
include only observations with CRF <10% in our sample. 
Uncertainty in the SCIAMACF1Y absolute calibration rela- 
tive to that of MODIS cause this check to allow a significant 
fraction of the data into our sample. Application of more or 
less stringent limits on cloud contamination within a moder- 
ate range did not substantially alter the derived spatial pat- 
terns of F. 

We calculated the fraction of a field-of-view for which the 
satellite does not sense the surface (/ e ff) for various cloud 
conditions. We can express / e ff as a linear function of the 
geometrical cloud cover / g , i.e., f e s = A x / g , where A is 
a parameter that depends upon surface and cloud optical 
properties. For a surface reflectivity of 0.2, A = ~0.5 for a 
cloud optical thickness (r) of 10 (using the Cl cloud model) 
and A = 0.75 for r = 20 at 760 nm (results will be similar at 
866 nm). Increasing surface reflectivity to 0.4 gives A = ~0.4 
for r = 10 and A = 0.75 for r = 20. This implies that in par- 
tial cloud cover with typical values of r (~10), a significant 
fraction of the surface is still seen by satellites. For example, 
in 40 % cloud cover with r < 10 and for surface reflectivities 
typical of vegetated land, 80 % or more of the surface signal 
is still seen. 

Owing to uncertainty in the absolute calibration of 
GOSAT, here we use the geometrical cloud fraction de- 
rived from the GOSAT CAI. We removed observations with 
/ g > 40 %. Derived spatial and temporal variations of ad- 
ditive signals did not depend substantially on the exact 
values of the thresholds used here for either GOSAT or 
SCIAMACHY. 

5.5 Quality control 


where f g is the geometric cloud fraction, and 7 0 bs and 7 c id 
are the observed and bare-cloud (no surface reflection) ra- 
diances, respectively. Flere we use the so-called mixed 
Lambertian-Equivalent Reflectivity model to estimate the 
CRF. In this model, 7 c id is set equal to / c id80, defined as a 
Lambertian surface with an equivalent reflectivity of 0.8, and 
/ g is set equal to an effective cloud fraction / e ff that is de- 
rived using the relationship 

fobs — /eff fcldSO T" (1 /eff) F\x- (3) 

(e.g., see Stammes et al., 2008, for an overview of the ap- 
proach that has been used to retrieve the cloud optical cen- 
troid pressure at NIR through UV wavelengths). In short, 
the selection of 0.8 as the cloud reflectivity produces an ef- 
fective cloud fraction that is generally less than the geometric 
cloud fraction and represents the fraction of the footprint that 
does not see the ground. Absorption and/or scattering from 
beneath a thin cloud is accounted for by the second term in 
Eq. (3). 

Here, we neglect atmospheric scattering as it is very small 
at 866 nm. To compute 7 c i r , we use the black-sky 16- 
day gridded filled-land surface albedo product from Aqua 


To create gridded data sets, we tested different filtering 
schemes for quality control. For SCIAMACHY, we included 
all data passing gross quality assurance checks on the filling- 
in, radiance residuals, and wavelength shifts. These checks 
remove very few observations. However, we have found evi- 
dence that the memory effect, described in Appendix A, pro- 
duces anomalous filling-in observed near coastlines where 
radiances can have sharp gradients. To eliminate the af- 
fected observations, we apply a filter that first compares a 
field-of-view continuum radiance with the average of the 
15 previously-observed fields-of-view continuum radiances. 
If the previously-observed average differs from the radiance 
by more than 15%, that field-of-view is discarded from the 
sample. This filter eliminated a significant amount of data 
near coastlines and thus results in somewhat noisier averages 
in those areas. For GOSAT, we removed outliers (>2 a) sep- 
arately for P and S polarization data as well as for the aver- 
aged P and S data. 


www.atmos-meas-tech.net/5/809/2012/ 


Atmos. Meas. Tech., 5, 809-829, 2012 


818 


J. Joiner et al.: Filling-in of near-IR solar lines 


5.6 Sealed fluorescence 

When comparing the derived additive signals with vegeta- 
tion indices, we use a quantity called “scaled-F” (Joiner et 
al., 2011), defined as the retrieved F divided by cos(SZA). 
This scaling roughly accounts for variations in F due to the 
incoming (clear-sky) PAR. 

5.7 Uncertainties 

Uncertainties for GOSAT retrievals were calculated by Joiner 
et al. (2011) using standard error propagation and assuming 
only random errors (i.e., instrumental noise). Here, we sim- 
ilarly compute errors for SCIAMACHY retrievals using the 
same assumptions. As discussed in Joiner et al. (2011), er- 
rors in scaled-F are proportional to R/SNR, and the relative 
error in terms of fraction of radiance is inversely proportional 
to SNR. For SNR= 1000, which is consistent with radiance 
residuals, we obtain a relative error of 0.3%. For R=Q2 
and 0.4, we obtained scaled-F errors of 0.37 and 0.75, respec- 
tively. Standard deviations in monthly scaled-F retrievals 
within gridboxes of 0.5° x 0.5° are roughly consistent with 
these expected errors as shown in Appendix C. 

However, errors from unaccounted-for systematic instru- 
mental errors can be quite significant. For example, the 
absolute GOSAT F values in this work are smaller (by 
approximately a factor of 2) than those derived by Joiner 
et al. (2011), owing to the new corrections applied. Be- 
cause unaccounted-for systematic errors in SCIAMACHY 
data could be significant and their magnitudes are unknown, 
it is difficult to accurately quantify the uncertainties in the de- 
rived scaled-F. Many SCIAMACHY observations are avail- 
able for averaging as compared with GOSAT which helps re- 
duce random errors. However, systematic errors are expected 
to remain. In addition to residual errors left over from imper- 
fect corrections, the South Atlantic Anomaly (SAA) has been 
known to affect measurements over South America. How- 
ever, we find no clear evidence of the effects of the SAA on 
our results. 

Here, we provide standard errors for gridbox averages to 
indicate uncertainty in the mean values. The number of ob- 
servations in a 0.5° x 0.5° gridbox in a month varies from 
0 to ~12. For the 8yr climatological monthly averages, 
the number of points in a gridbox thus ranges from 0 to 
about 100. Therefore, for gridboxes with relatively large 
numbers of observations, the standard error for the clima- 
tological monthly mean can be approximately an order of 
magnitude less than the standard deviation. Sample global 
maps of these statistical parameters are shown and further 
discussed for January and July in Appendix C. We note that 
these values may still underestimate errors in the derived ab- 
solute values of F and scaled-F and that theoretical error cal- 
culations based on the assumption of random, Gaussian er- 
rors may be considered only as a lower limit. 


6 Results and discussion 

6.1 Spatial variations in scaled-F 

Figure 8 shows retrieved gridded monthly mean scaled-F for 
July and December 2009 derived from GOSAT at 770 nm and 
758 nm and from SCIAMACHY at 866 nm. The values from 
SCIAMACHY at 866 nm are significantly smaller than those 
at 770 nm, as would be expected if the signals originate from 
chlorophyll a fluorescence associated with the declining 
emission tail throughout the NIR region. GOSAT results are 
comparable in magnitude and show similar spatial/temporal 
characteristics as compared with those of Frankenberg et al. 
(2011a) and Guanter et al. (2007). 

Similar seasonal variation in scaled-F is seen by both sen- 
sors as well as EVI, namely, higher Northern Hemisphere 
values in July versus higher values in the Southern Hemi- 
sphere in December. Similar spatial patterns are shown in 
scaled-F by both SCIAMACHY and GOSAT, such as low 
values over central Brazil in July and over the northeast of 
Brazil in December. These variations in GOSAT and SCIA- 
MACHY scaled-F are consistent with a vegetation source 
such as fluorescence. Like the EVI, a retrieved fluorescence 
signal (scaled by incoming PAR) is sensitive to the amount 
of green biomass contained within the sensor field-of-view 
or fractional amount of intercepted PAR. 

Subtle differences in spatial variations are shown be- 
tween GOSAT- and SCIAMACHY-derived scaled-F. SCIA- 
MACHY shows somewhat lower values at high northern lat- 
itudes in the boreal summer as compared with tropical re- 
gions. There could be several reasons for this if the pri- 
mary source of the SCIAMACHY signal is vegetation flu- 
orescence as is presumably the case for the GOSAT obser- 
vations. Firstly, there may be remaining errors due to instru- 
mental effects. Secondly, SCIAMACHY and GOSAT obser- 
vations are made at different local times; SCIAMACHY ob- 
servations are made near 10:00 LT while those of GOSAT are 
made about 2.5 h later, just after noon. The shapes of the flu- 
orescence spectra may differ for different types of vegetation 
as well as different magnitudes and spectral shapes of the in- 
coming solar light. For example, more UV light is absorbed 
in relatively clear skies at the surface in the tropics where 
total column ozone and SZAs are low. Plants at these lati- 
tudes have more photo-protective mechanisms in place that 
may produce different fluorescent spectral responses. Com- 
pounds within the vegetation besides chlorophyll a may al- 
ter the spectral dependence of fluorescence. Finally, canopy 
structure may influence the fluorescence spectra. Multiple 
scattering within a deep canopy may enhance fluorescence, 
producing a different spectral response than a sparse canopy 
(Rosema et al., 1991). 

Figure 9 shows monthly climatologies of scaled-F and 
pcos(SZA) (proportional to radiance) derived from over 8 yr 
of SCIAMACHY data (January 2003-December 2011, ev- 
ery other month, gridded at 0.5° x 0.5° resolution) and EVI 
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Fig. 8. Derived monthly averages for July (left panels) and December (right panels) 2009; first(second) row: fluorescence scaled by cos(SZA) 
(scaled-F) (unitless) from GOSAT 758 (770) nm; Third row: scaled-F from SCIAMACHY 866 nm; bottom row: Aqua MODIS enhanced 
vegetation index (EVI) 


from Aqua MODIS (August 2002-July 201 1 at the same spa- 
tial resolution) for comparison. An animated version of the 
SCIAMACHY scaled-F climatology is available in the sup- 
plemental material. The general spatial and temporal patterns 
from SCIAMACHY scaled-F are similar to EVI and consis- 
tent with fluorescence related to vegetation. Note that we 
show scaled-F here for solar zenith angles up to 80°. In the 
winter hemisphere at the highest latitudes shown, scaled-F 
originating from a terrestrial source is likely to be overesti- 
mated owing to atmospheric RRS that is not fully accounted 
for. 

Some interesting differences between scaled-F and EVI 
are shown here as well as in Fig. 8. For example, we 


see moderate values of EVI over parts of the northeast and 
southeast US, Europe, and southeast Asia in the late autumn 
through early spring while scaled-F is close to zero. The val- 
ues of scaled-F are more distinctly lower over parts of east 
and northeast Brazil from May through November as com- 
pared with the EVI. These differences are further explored in 
the next subsection. 

Scaled-F is shown to be loosely correlated with radiance 
in vegetated areas. This would be expected if the additive 
signal at the Ca II line has a vegetation source. However, 
there are also areas of high radiance with low F such as over 
the Sahara and parts of the Saudi Arabian peninsula. This 
demonstrates that the additive signal is not entirely related 
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Fig. 9. Monthly mean scaled-F (left panel) and p cos(SZA), proportional to radiance (right), derived from 8 + yr of SCIAMACHY data 
gridded at 0.5° x 0.5° resolution and gridded EVI from 9yr of Aqua MODIS at the same spatial resolution (middle panel) for every other 
month January through November (top to bottom panels). 










to radiance as might be the case if the signal was entirely 
due to instrumental artifacts. This is also consistent with a 
vegetation source for the 866 nm additive signals. Finally, we 
note that the retrieved value of scaled-F for SCIAMACFIY 
depends upon satellite viewing geometry (see Appendix C 
for details). 

6.2 Seasonal variations in SCIAMACHY and GOSAT 
scaled-F and EVI 

In Fig. 11, we examine the climatological seasonal cycle of 
scaled-F from both SCIAMACHY and GOSAT as well as 
EVI over several geographic regions shown in Fig. 10. For 
clarity, GOSAT 770 (758) mn values are reduced by constant 
factors of 2.25 (4.5). We show examples for several different 


biome types, using the same vertical range in all panels to 
highlight differences between scaled-F and EVI for the dif- 
ferent regions and vegetation types. A frill 8 yr time series for 
three of the regions with substantial seasonal cycles is shown 
in Appendix C. 

Over the deciduous broadleaf forests of the northeast US 
(Fig. 11, panel 1), SCIAMACHY scaled-F shows seasonal 
variations similar to those of GOSAT. SCIAMACHY val- 
ues in this region are proportionally smaller than GOSAT 
as compared with the other panels. Like scaled-F, the EVI 
shows a sharp rise starting around May and a decline start- 
ing in September. For evergreen needleleaf (pine) forests in 
the southeast US (panel 2), seasonal variations are similar 
for GOSAT scaled-F and EVI. SCIAMACHY scaled-F also 
shows a later springtime rise and and earlier autumn decline. 
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Fig. 10. Regions used for calculating seasonal cycles shown in Fig. 11. 


1) NE US: deciduous broadleaf 


2) SE US: evergreen needleleaf 



3) C.Amazon: evergreen broadleaf 
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Fig. 11. Seasonal cycles derived for the regions shown in the Fig. 10 map, using 0.5° x 0.5° gridbox monthly means within the indicated 
region and with the dominant vegetation type as listed: Black solid line: scaled-F derived from SCIAMACHY (left axis, unitless) and scaled 
GOSAT (red: 770 nm scaled by 2.25; blue 758 nm scaled by 4.5); vertical bars indicate error of the mean; Green dotted line: MODIS 
Enhanced Vegetation Index (EVI, right axis, unitless). The numbers indicated are the number of monthly mean 0.5° x 0.5° gridboxes used 
in each plot. 


For these and many other moderate to high latitude regions, 
values of the EVI have minima between 0.2 and 0.5 during 
winter months, while scaled-F is near zero. 

In the heavily forested areas of the central Amazon 
(panel 3), both the EVI and scaled-F remain relatively high 
year round. This contrasts sharply with shrubland over north- 
east Brazil (panel 4) at similar latitudes. Here, scaled-F and 
EVI are low in July through September following the green- 
up when scaled-F and EVI are high in January through April. 

The savannas of central Africa (panel 5) have a distinct and 
similar seasonal cycle in scaled-F and EVI. Scaled-F from 
SCIAMACHY shows an earlier peak than that of GOSAT 
over croplands of India (panel 6). GOSAT scaled-F also 
peaks somewhat earlier than EVI. The differences between 


scaled-F and EVI shown here remain unexplained; it is hoped 
that tli tu re studies will focus on this area to discern whether 
these differences stem from algorithmic/instrumental arti- 
facts or a geophysical effect. 

Although we have shown some differences between 
EVI and scaled-F, the similarity of the SCIAMACHY and 
GOSAT scaled-F temporal and spatial variations and the 
comparisons with EVI appear to be consistent with a major 
source of the SCIAMACHY scaled-F signals linked to veg- 
etationi fluorescence. The spatial and temporal patterns are 
not consistent with fire as a dominant source, though we note 
that fires may contribute to the observed signals. 
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6.3 Filling-in not consistent with vegetation 
fluorescence 

We note filling-in (retrieved as F) at 866 nm over parts of 
the Sahara desert (most apparent in May-September) and the 
Saudi Arabian peninsula where vegetation is sparse and a sig- 
nificant signal from chlorophyll fluorescence is not expected. 
We also note some negative filling-in values over these re- 
gions that have high radiances values. The memory effect of 
SCIAMACHY may contribute to these features as radiances 
at 866 nm show sometimes large gradients over these areas 
and our filtering scheme may not have completed removed 
all affected pixels. The SCIAMACHY memory effect may 
also affect the reference spectra, particularly for high values 
of radiance that correspond to the bright cloudy data over 
ocean. In addition, residual filling-in over the brightest areas 
may not have been well characterized by the cloudy ocean 
spectra. Cloudy spectra with high radiance values occur pri- 
marily in deep and frontal convection where cloud pressures 
are low. This may lead to an underestimation of filling-in 
from RRS when these spectra are used as a reference and 
thus over-estimation of the residual filling-in from terrestrial 
sources. 

Filling-in over barren regions may also be produced by lu- 
minescent minerals in soil and/or rock. Mineral detection 
was a primary goal of early measurements using the Fraun- 
hofer line discriminator method (e.g., Hemphill and Vick- 
ers, 1966; Watson, 1981; Stacy et al., 1984; Hemphill et al., 
1988). There are many types of minerals with impurities 
that luminesce when excited by UV and visible light (Gaft 
et al., 2005). For example, samples of fluorite (CaF 2 ) con- 
taining Nd 3+ (Neodymium) centers show laser-induced lu- 
minescence with peaks at 866 and 795 nm (Gaft et al., 2005). 
Similarly, Dy 3+ (Dysprosium) centers give narrow peaks at 
765 and 778 nm. 


7 Conclusions and ongoing work 

Our simulations indicate that terrestrial fluorescence filling- 
in of the 866 nm Ca II line can be detected using hyperspec- 
tral instruments (spectral resolutions of the order of tenths 
of a nm) such as SCIAMACHY if the fluorescence at this 
wavelength is of the order of 0. 1-0.2 mWm -2 nm -1 sr -1 . 
Laboratory measurements of vegetation fluorescence excited 
by broad-spectrum light indicate that such values are possi- 
ble. After corrections for instrumental artifacts, we retrieved 
an additive signal over land at 866 nm with SCIAMACHY. 
The magnitude of the derived additive signal at 866 nm is 
similar to that of our solar-simulator laboratory measure- 
ments. The spatial and temporal patterns of the detected 
additive signals at 866 nm are consistent with a vegetation 
source; they are similar to those of EVI and those derived 
from additive signals at 770 and 758 nm where fluorescence 
from chlorophyll a in vegetation is stronger and expected 


to be the primary source of the signals. Further laboratory 
and ground-based studies are needed to confirm that the ad- 
ditive signal retrieved at 866 nm is indeed primarily due to 
vegetation fluorescence. Measurements of fluorescence and 
reflectance spectra for a variety of vegetation types, canopy 
structures, illumination, and environmental conditions will 
help to improve the interpretation of satellite-derived addi- 
tive signals shown here. Additional sources of the additive 
signals include unaccounted-for rotational-Raman scattering 
in the atmosphere, fluorescence from minerals in rock and 
soil, and emission from flaming fires and volcanoes. 

Owing to the large number of observations collected over 
more than 8 yr of operation, SCIAMACHY data can be av- 
eraged to provide climatological spatial patterns of F at a 
higher resolution than is possible with the current set of 
GOSAT observations, even though the intrinsic spatial res- 
olution of SCIAMACHY is lower than that of GOSAT. Al- 
though the derived spatial patterns and seasonal cycles in 
scaled-F are similar to those of the MODIS EVI, impor- 
tant differences between scaled-F and EVI are noted. For 
example, a more distinct shut-off of activity in winter and 
dry months is shown in scaled-F in some areas as compared 
with EVI, which could improve estimates of the length of the 
physiologically-active growing season. 

More work needs to be done to further refine F retrievals 
from satellite instruments that were not designed to mea- 
sure small additive signals. Owing to the extremely small 
F signal derived from instruments that have relatively large 
systematic errors, our reported absolute values and tempo- 
ral/spatial variations may still contain residual errors. How- 
ever, the approach we developed appears to be capable of 
detecting and correcting for complex, time-varying instru- 
mental artifacts. 

In addition, more research is needed to fully exploit these 
data within global carbon models and other applications. For 
example, the influences of canopy structure and other param- 
eters, such as chlorophyll content, need to be further inves- 
tigated in order to understand their effect on large field-of- 
view satellite data (e.g., Rosema et al., 1991; Louis et al., 
2005; Meroni and Colombo, 2006; Middleton et al., 2008; 
van der Tol et al., 2009; Zarco-Tejada et al., 2009; Rascher 
et al., 2009; Damm et al., 2010). It has been suggested that 
when vegetation is not under stress and the variations in flu- 
orescence yield are small, fluorescence flux from the 740 nm 
emission peak, that is affected by reabsorption of emittance 
from the red peak, is driven primarily by PAR absorption 
(Moya et al., 2004). Another important factor is that the ra- 
tio of fluorescence from the 740 and 680 nm peak emission 
features expresses information about photosynthetic function 
and chlorophyll content (e.g., Lichtenthaler , 1987; Gitelson 
et al., 1999; Saito et al., 1998; Louis et al., 2005) and function 
(e.g., Middleton et al., 1996). Therefore, additional measure- 
ments of fluorescence at other wavelengths will be useful for 
retrieving more complete spectral emission curves that en- 
able improved relationships to vegetation function. In order 
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to more fully exploit satellite-derived fluorescence measure- 
ments for deriving relationships with key carbon-related pa- 
rameters such as gross primary production (GPP) or light-use 
efficiency (LUE), fluorescence retrievals should be combined 
with other observable quantities such as reflectance-based 
vegetation indices including the narrow band index based on 
two green bands, the PRI (Daumard et ah, 2010; Middleton 
et ah, 2011). 

The upcoming launch of the Orbiting Carbon 
Observatory-2 (OCO-2) (Crisp et ah, 2004) will facil- 
itate the use of solar Fraunhofer lines for fluorescence 
retrievals at wavelengths near the O 2 -A band (Frankenberg 
et ah, 2011a) and provide more opportunities to improve 
fluorescence measurements in the NIR region. The FLuo- 
rescence EXplorer (FLEX) (Rascher, 2007; European Space 
Agency, 2008), an ESA Explorer 8 Mission selected for 
Phase A/Bl in early 2011, plans to utilize the O 2 -A and B 
bands for chlorophyll fluorescence retrievals (Guanter et ah, 
2010) and other spectral information across the visible-NIR 
spectrum. FLEX would also measure fluorescence at a 
higher spatial resolution than current satellite sensors that 
were not designed or optimized for fluorescence retrievals. 

Appendix A 

Instrumental effects in SCIAMACHY data that may 
produce false filling-in 

Al Memory effect 

In channels 1-5, it was found that the signal deviates from a 
linear response. The deviation depends upon the signal level 
of the previous readout. It also depends on the signal level 
including analog offset and dark current. Since there is a de- 
pendence on the previous signal for a large range of detector 
fillings, artificial spectral features can be introduced into the 
measurements. 

We found systematic biases in Fraunhofer line filling as 
a function of position within a block of nadir measurements. 
The first several positions (following the limb-viewing mode) 
have the largest biases. Because of repeating orbital patterns, 
if this effect is not accounted for, a cross-track striping will 
be present in the observed filling-in, even after averaging a 
month of data. 

A2 Dark current 

There are two components that contribute to the dark signal 
in channels 1-5: the analog offset and the leakage current. 
The former does not depend on integration time while the 
latter does. We found different biases in the filling-in for 
the two primary pixel exposure times (PETs) used for nadir 
measurements in channels 4 and 5 that may be a result of 
unaccounted-for dark signals. 


A3 Wavelength calibration and undersampling 

Due to the Doppler and other effects, there are small wave- 
length shifts in Earth spectra with respect to the measured 
solar spectra. We found that these wavelength shifts vary 
systematically with latitude. Wavelength shifts between ob- 
served and reference spectra must be accurately accounted 
for when deriving a relatively small filling-in signal. 

A4 Stray light 

Stray light may be classified as two types: spectral and spa- 
tial. Both of these can distort observed spectra. Spatial 
stray light mainly affects the limb measurements of SCIA- 
MACHY, but may also produce errors in a fluorescence re- 
trieval in the nadir mode. Spectral stray light may be caused 
by diffuse reflection that adds signal to all detector pixels in 
a given channel, thus mimicking fluorescence. The largest 
expected uniform stray light fraction in SCIAMACHY is 
0.14% (Lichtenberg et al., 2006). 

So-called “ghost” stray light is caused by reflection from 
one wavelength to another. It does not add signal to all pixels 
equally, but can distort the spectrum. The correction method 
is accurate to around 25 %, leaving at most 1 % stray light in 
the spectrum. 

A5 Calibration 

Various issues with SCIAMACHY’s radiometric calibration 
are discussed by Lichtenberg et al. (2006). These may af- 
fect the absolute values of derived fluorescence. In addition 
to radiometry, the applied calibration scheme may introduce 
spectral features and thus produce false filling-in. 

Appendix B 

Temporal and spatial variations in GOSAT cloudy ocean 
spectra 

Here, we show temporal and latitudinal variations in the 
filling-in over cloudy ocean scenes caused primlarily by the 
zero-level offset problem in GOSAT data. Figure B1 shows 
filling-in as a function of normalized radiance binned at in- 
tervals of 0.05. We compute the filling-in for diagnostic pur- 
poses only; it is more straight-forward to interpret and com- 
pare with previous works than reference spectra. The overall 
dependence of filling-in on radiance is similar to the depen- 
dence derived using data over Antarctica by Frankenberg et 
al. (2011b). However, here we show that the filling-in at low 
to moderate radiances for both polarizations varies with both 
time and latitude. 
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Fig. Bl. Diagnostic filling-in derived over cloudy ocean pixels versus pcos(SZA) derived for GOSAT 758 (770) nm in left (right) panels 
for P(S) polarization in top (bottom) panels for July and December 2009 averaged over normalized radiance bins of 0.05 and separately for 
Northern (Southern) Hemispheres shown by solid (dashed) lines. 


Appendix C 

Parameters related to uncertainty in climatological 
monthly-mean scaled-F from SCIAMACFIY 

Figure C 1 shows examples of three mapped statistical param- 
eters (number of observations averaged, standard deviation, 
and standard error) for SCIAMACFIY gridded data in Jan- 
uary and July. All parameters correspond to monthly data 
obtained over the years 2003-2011. These data show that 
the standard deviations are roughly consistent with those that 
would be expected from instrumental noise. Additional dis- 
persion in scaled-F may result from variations in aerosol and 
cloud cover, actual variations in scaled-F owing to environ- 
mental stress on vegetation, and residual uncorrected system- 
atic errors. For most areas, the standard errors are substan- 
tially smaller owing to averaging of a number of observations 
taken over the 8 yr time period. 

For reference, Fig. C2 shows the entire time series of 
monthly-mean scaled-F for 3 of the boxes in Fig. 10 where 
there is a significant seasonal cycle. Inter-annual variabil- 
ity is seen in these data. Evidence of long-term trends due 
to either geophysical changes or calibration drift is not ob- 
vious. Future studies will focus on determining whether the 
inter-annual variations shown can be related to geophysical 
parameters such as environmental stress on vegetation. Be- 
cause the systematic errors that we attempt to correct for can 
represent a large fraction of the observed additive signal, we 
do not rule out the possibility that at least some of the inter- 
annual variations are due to uncorrected systematic errors. 
That is why we focus on climatological averages in this work. 


Finally, we note that there is a dependence of the SCIA- 
MACHY retrieved additive signal on viewing geometry. Fig- 
ure C3 shows scaled-F retrieved on the east and west sides 
of the SCIAMACHY swath, the difference between the two, 
and the east to west side swath difference of radiance for 
the July climatology. The regions with large cross-track 
asymmetry in scaled-F vary with season (not shown). The 
scaled-F angular dependence shows some similarities to that 
of radiance, but also shows some differences. Guanter et 
al. (2012) also report on viewing-angle dependence of re- 
trieved fluorescence from GOSAT. We did not obtain a sta- 
tistically robust dependence on viewing angle for GOSAT 
scaled-F retrievals. The observed cross-track asymmetry 
in SCIAMACHY scaled-F may be related to differences in 
the amounts of absorbed PAR and fluorescence escaping the 
canopy for different viewing geometries, related to leaf ori- 
entation and canopy effects. The results shown here are pre- 
liminary; further examination of these effects will be the 
topic of future studies. 

Supplementary material related to this article is 
available online at: http://www.atmos-meas-tech.net/5/ 
809/2012/amt-5-809-2012-supplement.zip. 
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Fig. Cl. Number of SCIAMACHY observations passing our quality control filters (top panels), the gridbox standard deviations (middle 
panels), and standard errors (bottom panels) for the January climatological average over years 2003-201 1 (left panels) and the same for July 
(right panels). 
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Fig. C2. Monthly mean scaled-F as derived from SCIAMACHY versus time over the period January 2003-December 2011 for 3 of the 
boxes shown in Fig. 10 as indicated. Standard deviations are shown with red vertical bars and standard errors are shown with thick black 
vertical bars (sometimes too small to be noticeable). 
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Fig. C3. Monthly mean scaled-F as derived from SCIAMACHY for observations on the left side of the swath (top left panel) and right side 
(bottom left panel) for July averaged over the years 2003-2011. The difference is shown in the top right along with the similarly derived 
difference in radiance (bottom right panel). 
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